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Charmonium production and absorption in heavy ion collisions is studied 
with the Ultrarelativisitic Quantum Molecular Dynamics model. We compare 
the scenario of universal and time independent color-octet dissociation cross 



t""- , sections with one of distinct color-singlet J/V', ip' and Xc states, evolving from 

small, color transparent configurations to their asymptotic sizes. The mea- 

fSJ ' sured J/tp production cross sections in pA and AB collisions at SPS energies 

^^ ■ are consistent with both — purely hadronic — scenarios. The predicted ra- 

Qv^ ■ pidity dependence of J/ip suppression can be used to discriminate between the 

(— I , two experimentally. The importance of interactions with secondary hadrons 

Qh! and the applicability of thermal reaction kinetics to J/ip absorption are in- 

Qh! vestigated. We discuss the effect of nuclear stopping and the role of leading 

hadrons. The dependence of the ip' /J/tp ratio on the model assumptions and 

the possible influence of refeeding processes is also studied. 
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I. INTRODUCTION 



Experimental data from the CERN SPS [|I]-|^ clearly show that the yield of J/ ^-mesons 
in pA and AB collisions is suppressed when compared to hard QCD production. A sub- 
stantial portion of these data can be described in simple models of nuclear absorption with 
certain assumptions about the dissociation cross sections, e.g. [^p!^. Recent, preliminary 
data on J/ip production in Pb+Pb collisions [0-^ have generated an intense and ongoing 
controversy regarding the necessity of a novel suppression mechanism, e.g. a quark-gluon 



plasma [|1J], to explain the data or whether the suppression can be due to hadronic interac- 
tions. (For theoretical reviews, see Refs. |TB|-P^].) Charmonium production and absorption 
has been studied first by means of semianalytical models of heavy ion collisions and/or ther- 
mal rate calculations |^-|T^,|TP|-^ as well as, more recently, within microscopic hadronic 
transport models p5|-p0|. Some current models suggest purely hadronic explanations of 
the suppression while others claim that this is not possible. Here we use different model 
scenarios to identify crucial theoretical and experimental conditions which have to be met 
before definite conclusions can be drawn. 

In simple analytical models of J/ip suppression, based on a semiclassical Glauber ap- 
proach, nuclear absorption is calculated by convoluting nuclear density profiles with an 
assumed cross section for dissociation by nucleons. The model parameters can generally be 



tuned to fit all pA data as well as all the AB data except Pb+Pb 0-0 . Absorption by pro- 



duced hadrons, comovers, is typically also included in these models [[7|,P|-[r^. The comover 
dynamics is treated in a simple, schematic fashion, generally assuming that the comover 
density scales with the number of nucleon participants with the exception of Ref. [l^. It is 



not straightforward to distinguish between these results because of the variation of model 
parameters such as the comover dissociation cross sections, nuclear density profiles, forma- 
tion and freeze-out times, as well as including feed-down to the J/ip from more massive 
charmonium states. However, except for Ref. |l^, all these studies claim that hadronic 



interactions alone cannot account for the J/ifj suppression in all pA and AB data. 

Implicitly, geometrical models of nuclear absorption such as Glauber scattering assume 
that the nuclei traverse each other without distortion or energy loss. This seems to be appro- 
priate for the description of the hard Drell-Yan |]3T|j3^ or cc [Q production processes and 
some authors also treat charmonium dissociation perturbatively [Q. However, the dissoci- 
ation of charmonium need not necessarily be associated with a large momentum transfer. 
Therefore, nuclear absorption may result from soft processes which have been shown to 
dominate the charmonium dissociation at SPS energies [^]. If nuclear absorption is a soft 



process, the nucleons may scatter and lose energy before interacting with a charmonium 
state. Microscopic transport models, which approximate the soft hadron dynamics of a 
heavy ion collision as a cascade of binary scatterings, include this nontrivial baryon dynam- 
ics and are therefore adequate tools for the simulation of charmonium absorption. Moreover, 
they treat the complete space-time evolution of secondary hadrons in a consistent fashion. 
However, the treatment of hard production processes as gg — ^ cc or qq -^ yU+/i~ is not 
obvious in the microscopic simulation of a soft hadronic cascade. The Drell-Yan production 



cross section in pA interaction is known to scale linearly with A ||31|J3^ : UpA = ctqA . Because 
of the apparent charmonium absorption discussed below there is no direct experimental 
evidence for linear scaling in charmonium production, or the hard nature of this process, 



respectively. However, we assume hard production of charmonium since measurements of 



open charm production imi + A reactions show no significant deviation from A scahng [^ 
It should be mentioned, though, that the open charm data at low transverse momentum 
suffers from low statistics. We also note that, even though charmonium production does not 
scale as A, the exponent, 0.9 [^,Q, is closer to unity than for lighter hadrons at low pt 
and xf which grow approximately as A^^'^ [p^j . 

The near linear scaling of the cc production in pA reactions is a consequence of requir- 
ing a single perturbative parton-parton interaction. To simulate this microscopically, soft 
interactions of nucleons and the associated energy loss must be neglected. In Refs. P^ , pB| , 
cc production was implemented into the hadronic transport code as a possible subchannel 
of the total A^A^ cross section. This is certainly not consistent with the picture of pertur- 
bative cc production. In Ref . [p^-|30[] , on the other hand, the space-time production points 
are obtained from hard processes in each event before the hadronic background simulation. 
Charmonium dissociation is also treated differently in the various models. Some simula- 
tions, p6|j28| , |29[| , assume that charmonium precursor states can be immediately dissociated 
by nucleons with constant cross sections larger than 5 mb. In Ref. [^] the charmonium 

color transparent states while Ref. 



cross sections evolve from small, color transparent states while Ref. I^S] gives no definite 
formation time. Only Refs. pU| , |5(]| include feeddown from ip' and Xc mesons. Based on 
rather different model treatments of the charmonium dynamics and the interconnection of 
hard and soft processes, all these studies claim that conventional hadronic scenarios are 
consistent with the preliminary Pb+Pb data. 

In this work, we simulate three different scenarios of charmonium absorption within 
the same calculational framework and compare the results. We apply the Ultrarelativistic 



Quantum Molecular Dynamics model vl.l [^ for all calculations presented here. 

In the first scenario, charmonium states are assumed to be produced as color-octet states 
which can be dissociated by nucleons only ||ll[. In the second scenario, nuclear absorption 
is treated in the same way as the first. However, secondary hadrons may dissociate the fully 
formed J/i/^, ip' and Xc states after a finite formation time. Finally, we investigate a scenario 
in which the charmonia are produced as distinct color-singlet states and evolve from small, 



color transparent configurations to their asymptotic sizes [ p3| , p0|1 . 

In Sec. |l|, the model framework and the distinct scenarios will be discussed in detail. 
The results of the simulations are then presented and discussed in Sec. |ITT[ In Sec. [111A| , we 
examine the consistency of transverse energy spectra measured by different experiments and 
compare them with our calculated distributions. In Sec. [ill B| , some qualitative features of 
charmonium suppression and the space-time evolution of the system are discussed. Sec. |IIIC| 
compares the quantitative AB systematics of J/ip suppression within our scenarios to the 
data. After a short discussion of the basic differences between semianalytical and micro- 



scopic descriptions of the comover dynamics in Sec. |III D| , we turn to the ijj' / J /ijj ratio and 
its sensitivity to different model assumptions in Sec. |11IE| . In Sec. |111 h\ the dependence of 
J /ip suppression on the projectile energy is studied. The key role of rapidity in discriminat- 
ing between the different scenarios J /ip suppression is identified in Sec. |illG . Finally, we 
investigate the time dependence of J / ip-covaoYei collision energies and dissociation rates in 
Sec. [1IIH| . The applicability of thermal reaction kinetics is also discussed. In Sec. 0, we 
summarize our results and draw conclusions. 



II. THE MODEL 
A. The general framework 

We simulate the production and the absorption of charmonium states in heavy ion coUi- 



sions within the Ultrarelativistic Quantum Molecular Dynamics model ||5^. This model has 
been shown to successfully describe final rapidity distributions of negatively charged hadrons 
and net-baryons in p + p and nucleus-nucleus collisions |^S[| . While the model prescription 
for the charmonium rescattering stage is varied in order to study different conceivable sce- 
narios, the initialization of the nuclei and the simulation of the hard cc production process 
are the same for all calculations presented in this work. The large number of parameters 



in the phenomenological microscopic model |^, especially in its collision term, have not 



been altered in the present study. The default values of these parameters are fit to data 
from elementary hadronic and p-|-A reactions over a range of energies. All concepts and 
parameters related to the treatment of heavy quarkonia in the model are specified in the 
following. 

We apply perturbative QCD to the production of charmonium states by simulating 
nucleus-nucleus collisions in the impulse approximation. This is done microscopically by al- 
lowing freely streaming projectile and target nucleons to scatter without deflection or energy 
loss. The nuclear dependence of the parton distribution functions is neglected. The result- 
ing space-time distribution of charmonium production points is inserted into the evolving 
hadronic environment calculated with UrQMD since the rare quarkonium production pro- 
cesses are only small perturbations on the heavy ion collision^. Our model is thus designed 
to account for partonic and hadronic aspects of the charmonium dynamics. 

In UrQMD, the initialization of the projectile and target nuclei is based on the Fermi- 
gas ansatz. The centroids of the nucleons are then randomly distributed within a sphere of 
radius 



RiA) = ro 



\'^ 



^1/3 _ 1 



1/3 



(1) 



38|| . The effective interaction radius is somewhat larger than R{A) since the individual 
nucleons have a finite size. In pA reactions, a smeared out surface region of AR =1.1 fm, 
extending beyond R{A), results from the geometrical A^A^ cross section, cr ~ 40 mb. The 

parameter tq = (4:;^— )^ = 1-14 fm since the nuclear matter ground state density po = 
0.16 fm~^ is used in UrQMD. The above described initialization renders a good agreement 



with measured charge density distributions in nuclei [p9|] . 

We use this initialization procedure for the simulation of hard processes via perturbative 
QCD and the soft hadronic cascade, respectively. Consequently, the average nuclear density 
profiles are the same. However, in each event, the initial positions of the nucleons are deter- 
mined independently for the two simulations after the impact parameter b is chosen. Thus, 



^ Since we consider only exclusive charmonium production, the QCD factorization theorem is 
inapplicable. 



we insert the cc states into the hadronic environment according to an averaged distribution 
of hard production processes rather than hnking the hard and soft simulations event by 
event. We therfore avoid strong and probably unphysical spatial correlations between hard 
and soft processes^. 

In the microscopic UrQMD transport model |3^ , many nucleons are excited to baryonic 



strings, i.e. their constituent quarks and diquarks are separated and propagate indepen- 
dently as the endpoints of color flux tubes. A significant fraction of the constituent quarks 
in the incident nucleons may thus appear as the core of leading mesons and are allowed to 
interact immediately after their production while leading baryons containing the remaining 
diquark interact with a reduced cross section in the UrQMD model. A sharp distinction 
between nuclear absorption and absorption by secondaries is not supported in this picture 
of nucleon break-up. 

B. Distinct scenarios 

We now introduce the scenarios for quarkonium suppression applied in our simulations. 

1. Scenario I 

This scenario follows Ref . |jll[ . Charmonium states are assumed to be produced as color- 
octet 'preresonance' states, \ccg), which exist until after their passage through the nuclei is 
completed. These preresonance states can be dissociated by nucleons with a constant and 
universal cross section. In the NN center-of-mass frame, the \ccg) states are produced at 
rest. Nuclear stopping and a change of the absorption cross section due to string excita- 
tion of incident nucleons is not considered. Then, in the UrQMD simulation, all but the 
charmonium-nucleon interactions have zero cross sections. J /ip dissociation by secondary 
hadrons is not taken into account, since QCD calculations with heavy quarks predict a 



strong damping close to threshold |Q. There is one free parameter, the |cc(7)-nucleon cross 
section. 



2. Scenario II 

According to this scenario, nuclear absorption is due to interactions of color-octet states. 
Thus, absorption by nucleons is calculated as in Scenario I, and the \ccg) dissociation cross 
section is the first free parameter. Thus, we follow the arguments of [jll| and the 'prereso- 
nance' hypothesis, respectively, as far as the nuclear absorption is concerned. 



^If the hard and soft simulations are linked event by event, the resulting absorption probabilities 
sensitively depend on technical details of the implementation, e.g. whether the position of a 
produced J /ip is chosen to be exactly the arithmetic mean of the positions of the two incident 
nucleons or randomly distributed between them. This seems to be unreasonable and is due to the 
naive geometrical interpretation of the interaction cross section in the semiclassical cascade. 



Then, in a full UrQMD simulation, the environment of secondary hadrons is calculated. 
After their finite formation time (on average, rj? ~ 1 fm/c) resulting from the string frag- 
mentation used in the UrQMD model, secondary hadrons may dissociate the J /'ip, ip' and 
Xc) which are assumed to be fully formed states by Tp- In order to clearly separate the nu- 
clear absorption of color-octet states from the delayed interactions of produced mesons with 
the asymptotic color-singlet states, we exclude the prompt interactions of leading secondary 
hadrons. The interaction of the leading quarks is effectively accounted for in the calculation 
of the nuclear absorption because the possible break-up of the incident nucleons due to a 
string excitation is ignored. 

To calculate charmonium dissociation by secondary hadrons, it is necessary to distinguish 
between the distinct charmonium states in the simulation. The charmonium states are 
distributed according to their assumed production probability times their decay probability 
to J/ip^s. Thus 40% of the final states are Xc's, 55% are J/ip^s, and 5% are ip's [|^. According 
to the degeneracy of the angular momenta, 1/3 of the Xc's are Xcio states and 2/3 are Xcii 



states. Their momenta are assigned according to the parametrization |41 



^d^ ~ (1 - ^F?-'' exp{-pT 2.08 GeV-i) . (2) 



Unlike Ref. [|r^, we assume that the comover dissciation cross sections are predominantly 



nonperturbative at SPS energies, as shown for the dissociation by nucleons pj]. The nucleon 
absorption cross sections, X{cc) + N, were derived in the framework of a nonrelativistic 



quarkonium model [^, obtaining a{J /^'N) = 3.6 mb, a{^jj'N) = 20 mb, cr{XcioN) = 6.8 mb, 
and a{xciiN) = 15.9 mb. These cross sections are not relevant for nuclear absorption in this 
scenario which is assumed to be exclusively due to interactions of the hypothetical color-octet 
states, as discussed above. The charmonium-meson cross sections {X{cc) + ii, X{cc) + p, 
etc.) are obtained by scaling the nucleon absorption cross sections by a common factor F 
so that a{X{cc)M) = F a{X{cc)N). In the additive quark model, F = 2/3. However, in 
this scenario it remains a second tunable parameter. Thus, additional suppression of the 
comover dissociation cross sections due to the relatively lower X{cc) + M collision energies 
can effectively be accounted for by adjusting F. All comover interactions above the respective 
dissociation thresholds are assumed to break up the charmonium state. 

Charmonium dissociation is thus calculated in two different simulation modes for nu- 
clear and comover absorption, respectively, assuming factorization of nucleon and comover 
absorption. Separate survival probabilities for dissociation by nucleons, Sn, and mesons, 
Sm, are computed. The total survival probability is then St = SnSm- There are similarities 
between this scenario and the two-stage cascade used in Ref. 123 . 



3. Scenario III 
Here, the charmonium states are produced as distinct color-singlet states. Their dis- 



sociation cross sections evolve linearly from small, color transparent configurations p3 
Charmonium absorption is considered as a nonperturbative process, thus nuclear stopping 
is included in the calculation of nuclear absorption. In particular, nucleons which are excited 
to strings in a prior NN collision have a reduced cross section since they exist only as leading 



diquarks during their formation time. Charmonium states can interact with incident nucle- 
ons, leading quarks and diquarks, and secondary hadrons (after their formation time). The 
production probabihties of the individual charmonium states and their momentum distribu- 
tions are chosen as in Scenario II. The rescattering cross sections for X{cc) + B, assuming 
B = N, are taken from Ref. |^ (see Scenario II). We now let F = 2/3 according to the 
additive quark model. All baryon and meson collisions above the respective dissociation 
thresholds are assumed to break up the charmonium state. The cross sections correspond to 
the geometrical transverse radii r^ = \/— of the charmonium states. We use a' to estimate 
the respective formation times Tp of the charmonium states by choosing Tp = rJj,/c. During 



these formation times the cross sections increase linearly with t [^, starting from ctq = 
at t = 0. Since all the parameters have been fixed either by prior model calculations p3[ or 



the dynamical model, UrQMD, itself, we do not adjust any parameters in this scenario. 

III. RESULTS AND DISCUSSION 

A. Et distributions 

The transverse energy distributions in different systems at SPS energies have been cal- 
culated with the UrQMD model |3^. Figure |1| shows the differential cross section of the 
transverse energy produced in Pb(160 GeV)+Pb collisions within the pseudorapidity range 
2.1 < ?7 < 3.4. The UrQMD result is compared with experimental data from the NA49 



Collaboration |^2|. The overall agreement appears to be good although the detailed shape 
of the spectrum at very central events, i.e. very high Ej-, is somewhat different. This dif- 
ference may be due to statistical limitations of the simulation since we have analyzed only 
1000 events. The disagreement between model and experiment at very low Et can be due to 
an impact parameter cut at 6 = 12.5 fm in the calculations which leaves out part of the very 
peripheral reactions. The comparison between the S(200 GeV)+Au data from the NA35 
Collaboration and the model calculation within the same acceptance, 2.1 < r^ < 3.4, shows 
an agreement of similar quality, see Fig. |l[ 

The hard collision spectrum of neutral Et for Pb(160 GeV)-(-Pb and S(200 GeV)+U 
reactions has also been calculated, based on the following concept: First, the number of 
hard nucleon-nucleon collisions, which is proportional to the yield of Drell-Yan muon pairs, 
is determined microscopically in the simulation of freely streaming nucleons, as described in 
Sec. [11 A|. The associated simulation of all soft processes with the full hadronic cascade, ini- 



tialized with the same impact parameter b, then renders the information about the produced 
transverse energy. 

Figure ^ (a) shows the number of Drell-Yan muon pairs (oc number of collisions in the 
simulation of hard processes) in Pb+Pb interactions as a function of the produced neutral 
transverse energy within 1.1 < 77 < 2.3. The normalization is arbitrary. Experimental 
data from NA50 [|^ have been put on the plot as well. However, the abscissa has been 
rescaled since the most recent data [ffl indicate that the absolute Et values of the previous 



publications on Pb+Pb collisions pj43|] have shifted. We are aware that the new analysis 
does not imply a simple overall linear rescaling of the old data points. However, in order to 
reasonably compare the gross features of the experimental dimuon Et spectrum with our 
model calculation, we have multiplied all iJ^-values of the data in the plot by 0.8. This 



value was estimated by comparing the Et — Ezdc contour plots of the Tsukuba '97 |^ and 
the Moriond '98 [|] proceedings. The UrQMD simulation seems to be in rough agreement 
with the resulting dN/dEx spectrum. The agreement between the model and experiment is 
not unexpected since the production of the NA49 transverse energy had been shown to be 
correctly described. The additional free streaming simulation is very simple and represents 
a well understood model of hard scattering processes in nucleus-nucleus collisions. 

The agreement between model and experiment (NA38) P],P[ seems to become poor in 
the case of the asymmetric S+U system, see Fig. |^ (b). The UrQMD calculation appears to 
significantly overestimate the neutral transverse energy within 1.7 < r] < 4.1 although the 
Ej- spectrum of the similar S+Au system agrees with the NA35 data, as was shown in Fig. |1|. 
Note that the use of a model for the simulation of hard collisions and the weighting of the 
Et distribution with the respective number of Drell-Yan dimuons affects only the shape of 
the Et spectrum as compared to a normal distribution of all events — central events are 
favored, since they correspond to a higher number of nucleon-nucleon collisions. Thus, even 
if the correlation between transverse energy and impact parameter was completely wrong 
in UrQMD, which seems unlikely given the agreement in Fig. |^, this would merely lead to 
an incorrect weighting of the respective Et values. In Fig. |^ (b), however, one sees that the 
maximally achievable neutral transverse energy differs by about 25%. 

In view of our simulation results, the NA38 neutral transverse energy spectrum 
of S(200 GeV)+U seems to be inconsistent with the transverse energy spectrum of 
S(200 GeV)+Au measured by NA35. On the other hand, the (rescaled) NA50 Pb+Pb 
spectrum appears to be consistent with the NA49 and the NA35 spectrum. Given the con- 
sistency of the microscopic simulation with three different experimental measurements we 
assign the remaining discrepancy in S-t-U to a possible problem in the determination of the 
NA38 Et scale. In order to compare our calculations of J /ip production as a function of Et 
with NA38 data, we will therefore apply a scaling factor of 1.25 to the Et values reported 
byNA38[|T],|l. 

B. Space-time picture of charmonium production and absorption 

In this section, we illustrate how charmonium production and absorption is incorporated 
into UrQMD. Figs. |^ (a)-(d) visualize the production and absorption points of cc states in 
the z-t plane for central Pb(160 GeV)+Pb collisions as calculated with UrQMD. For this 
qualitative analysis, time independent baryon and meson absorption cross sections have been 
used. Their absolute value is irrelevant for this discussion. The impact of a time dependence 
of the cross sections on the dissociation probabilities will be discussed later in this section. 

The straight lines indicate the light cone trajectories of the geometrical boundaries of 
the nuclei. Hard production processes occur only at times and positions where the nuclear 
densities of projectile and target actually overlap in the course of the evolution. Since the 
nuclei traverse each other without energy loss in the present simulation of hard processes, the 
overlap region is determined by the light cone bands in the plot. Figure |^ (a) shows the cc 
production points. They all lie in the nuclear overlap region. Figure ^ (b) shows the points 
of cc dissociation by 'hard nucleons' in the z-t plane, i.e. nuclear stopping is not included 
in the simulation of nuclear absorption. Therefore, nuclear absorption is constrained to the 
region covered by the light cone bands of the two nuclei, however not to the actual overlap 
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region of both bands. The charmonium states are produced at rest. Therefore, the spatial 
distribution of absorption points must coincide with the distribution of production points, 
i.e. it is restricted to its longitudinal range \z\ < 0.8 fm. The temporal distribution of cc 
dissociation, however, naturally exhibits some delay as compared to the production, since 
absorption does not take place instantaneously. 

Fig. 1^ (c) shows the points of cc dissociation by nucleons as in Fig. ^ (b), but the 
charmonium states are now produced according to their assumed momentum distribution 
given by Eq. (^. Therefore, the spatial coordinates of production and absorption do not 
coincide. However, due to their high mass, few charmonium states actually acquire velocities 
high enough to propagate for very long and still be within reach of the nuclei which did not 
suffer energy and velocity loss in this calculation. 

Fig. ^ (d) shows the points of cc dissociation by 'soft nucleons' in the z-t plane, i.e. 
nuclear stopping is included in this simulation of nuclear absorption. Nuclear absorption 
can therefore occur anywhere within the forward light cones of the two nuclei. Most of the 
nuclear dissociation processes, however, take place at very early times, t < 3 fm/c. 

Fig. ^ (a) shows the points of cc dissociation by mesons in the z-t plane for central 
Pb(160 GeV)+Pb collisions. Note that mesons containing a constituent quark of an incident 
nucleon are included in this plot. These leading hadrons can interact with other particles 
instantaneously, however, their cross section is reduced by a factor of two. Most of the 
absorption processes within the first 1.5 fm/c are due to interactions with leading mesons. 
To demonstrate this. Fig. | (b) shows the points of cc dissociation by mesons in the z-t plane, 
where only interactions of fully formed mesons are taken into account. Clearly, if leading 
mesons are considered as comovers, they dominate the comover absorption probability. Even 



a strong damping of the dissociation cross section close to threshold |Q would not justify 
neglecting these processes, since they occur well above threshold (see Sec. |llli:d| ). On the 
other hand, nuclear stopping reduces the nuclear absorption considerably. This is due to the 
fact that when a nucleon is excited into a string, it loses one of its constituent quarks which 
becomes part of a leading meson. Accordingly, the dissociation cross section of the baryon 
containing the remaining leading diquark is rescaled to 2/3 of the nuclear absorption cross 
section in the UrQMD simulation. 

C. AB systematic of J/ip suppression 

Here we discuss the systematics of J/tp production and suppression as a function of 
AB in our three scenarios. We compare our results to the data from several projectile- 
target combinations: p(450 GeV)+C, p(200 GeV)+Cu, p{200 GeV)+U, S(200 GeV)+U, and 
Pb(160 GeV)+Pb. The 450 GeV and 160 GeV simulations are rescaled to piab = 200 GeV 
with the parametrization of Ref. |45|, as done by NA50 |^J^. 



1. Scenario I 



Figure ^ shows a comparison between the semianalytical Glauber calculations of 
Kharzeev et al. [O] and UrQMD /Scenario I using the same absorption cross section. 



^diss _ ^ 2i mb. The experimental data are also shown [Q. The J/ip., ip' ^''^^ Xc are as- 
sumed to interact with nucleons while in \ccg) color-octet states. 

Charmonium dissociation by secondaries is assumed to be negligible. As in the geomet- 
rical picture, the cc states are produced at rest and the nuclei traverse each other without 
energy loss. Nucleons do not suffer a reduction of their cross section due to string exci- 
tations. The difference between the semianalytical model and the microscopic simulation 
according to Scenario I is 10% in S-(-U collisions (Sciauber = 0.49 vs. SurQMD = 0.45). This 
is due to the different nuclear density profiles: In Ref. |Tl| a Woods-Saxon shape is used, 
whereas UrQMD initializes hard spheres with an effectively smeared-out surface layer due 
to the finite geometrical A^A^ cross sections, see Sec. [11 A| . Note that the cross sections are 
derived by multiplying the respective survival probabilities with the experimental p+p cross 
section 0. In Ref. |11|, p+D reactions have been used as reference value for determining 
the nuclear absorption cross section. Surprisingly, the p+D point is above the p + p point, 
measured at the same energy. Note that Fig. ^ is not our best fit through the pA and S-l-U 



data. As in Ref. [|n[|, the present (analogous) calculation underpredicts the J/ip absorption 
in Pb+Pb, while reasonable agreement for all other systems is achieved. However, a value 
of 9 mb for the color-octet absorption cross section yields an acceptable simultaneous agree- 
ment with all data if the p + p cross section is left as a free parameter, as shown in Fig. |^. 
Note that our p + p value differs from the data by two standard deviations in this case. 

2. Scenario II 

Figure |^ shows a calculation of the J/ip systematics within Scenario II. Recall that in 
order to avoid an overlap between the interactions of the charmonia as color-octet and color- 
singlet states, leading mesons are not allowed to interact. The early stage of cc absorption 
is thus described solely by Glauber-like nuclear absorption using a universal color-octet dis- 
sociation cross section, a{X{cc)N) = 4.8 mb, cf. [0]. Interactions with secondary hadrons 
occur only after their formation time — on the average r^ ~ 1 fm/c — after the charmonium 
states have fully formed. For the comover absorption cross sections we take F = 0.5 with 
the nucleon absorption cross sections from Ref. [^. The calculated J/ip production cross 
sections for the given parameter values all lie within the experimental error bars. The mea- 
sured cross sections are therefore consistent with a purely hadronic picture of the dynamics 
as described by Scenario II. Note that similar model assumptions in a different hadronic 
cascade model with a somewhat different choice of parameters have lead to the same con- 



clusion ||2^. Semianalytical calculations, however, employing the same nuclear dissociation 
cross section [^ fit the S+U data with significantly smaller comover absorption cross sec- 
tions. In contrast to the microscopic calculation presented here, simultaneous agreement 



with the Pb+Pb data cannot be achieved in that analysis, as will be discussed in Sec. [IIID . 



3. Scenario III 

Figure ^ shows the UrQMD calculation assuming the charmonium states are produced 
as color-singlet states and evolving from small, color transparent configurations. Nuclear 
stopping and dissociation by leading diquarks and quarks is taken into account. One observes 
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a fair agreement of the simulation with the AB systematics of the data. The measured cross 
sections are therefore also consistent with a purely hadronic system. Note that no parameters 
were tuned to fit the data. 

Figure ^ shows the resulting Ex distributions for the ratio of J/ip to Drell-Yan production 
in Pb+Pb and S+U reactions. The gross features of the Et dependencies of the ratios are 
reasonably well described by the model calculations. In particular, the calculated result is 
smooth in both cases, without abrupt discontinuities, in agreement with new high statistics 
data comparing J/ip production to the minimum bias cross section [^. Apparent disconti- 
nuities in the preliminary data which compared the more abundant J/ip to limited statistics 
Drell-Yan production had been interpreted by some as a sudden onset of a deconfined phase 



We want to point out that possible agreement or disagreement of the simulation with 
the final data depends decisively on the determination of the correct E^ scale. When deal- 
ing with preliminary data, horizontal error bars indicating the systematic uncertainties in 
the Et scale of the experimental analysis would be helpful. The uncertainties of the E^ 
determination in our calculation are only statistical and lie within the bin width of the 
histograms. 

Note that the slight underestimation of the J/tp over Drell-Yan production at low E^ 
seen in Fig. ^ can be related to an overestimation of the Drell-Yan cross section at low Et, 
visible in Fig. ^ (a). Calculations in the context of the Dual Parton Model also show have 



feature fl3 



D. Semianalytical vs. microscopic descriptions 

It has been shown in the previous analysis that, within the microscopic hadronic trans- 
port model UrQMD, comover interactions can explain the observed J/ip suppression in S-l-U 
and Pb-|-Pb collisions. The question then naturally arises why semianalytical calculations of 
charmonium production and absorption, based on the Glauber approach and the wounded 
nucleon model, cannot also simultaneously describe these data [|TO|JT^ . Therefore, it is worth- 
while to examine some basic differences in the comover dynamics of the two approaches. 



First, transverse expansion is neglected in the semianalytical models [[r0|-p!3|] although it is 
a feature of the true reaction dynamics, as established experimentally |^ . However, the full 
system including the transverse dynamics ||50| is treated explicitly in microscopic models like 
UrQMD. It remains to be seen how the J/ip survival probabilities are quantitatively affected 
by the inclusion of a realistic transverse expansion in the semianalytic models. 

Secondly, comover scenarios based on the wounded nucleon model assume a linear scaling 
of the produced hadron rapidity densities with the number of participants, Ed^N/dp^ ~ 
-^part where a = 1. In contrast, a possible enhancement of the meson yield at midrapidity 
for nuclear collisions relative to N+N collisions has been observed [^. The h~ yield per 
participant nucleon pair, extrapolated to full phase space, is 4.0 ±0.2 for Pb+Pb, compared 
to 3.6 ± 0.2 for S+S [Q. Measurements of the neutral pion production in Pb+Pb collisions 
as a function of centrality show that the pion yield in the interval 0.5 < pt < 3 GeV/c 
scales like iVirt @' Thus, experimental observations suggest deviations from the simple 
assumptions of the wounded nucleon model at SPS energies. This is not surprising since 
measured particle momentum correlation patterns are in striking contrast to simple p + p 
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and p + A extrapolations due to additional meson-meson and meson-baryon interactions in 
nucleus-nucleus collisions |54 . 



In UrQMD simulations of central Pb+Pb relative to central S+U collisions the meson 
yield scales like N^f^^. This nonlinear scaling with the number of participants naturally 
explains at least a part of the "additional" suppression in the heavier system. Note that 
semianalytical studies where the comover densities have been extracted from the Dual Parton 
Model, which also exhibits a nonlinear scaling, produce simultaneous agreement with the 
S+U and Pb+Pb data [|13[. 



Finally, semianalytic models generally neglect the considerable rapidity asymmetry of 



hadron production in asymmetric nuclear collisions ||5^ while UrQMD simulations also 
show this feature of the measured particle distributions. Even if the yields of produced 
hadrons over Air scaled exactly with the number of participants, the relative comover densi- 
ties extracted from the wounded nucleon model would be incorrect in the finite experimental 
acceptance cut of NA50 , < y < 1. According to the UrQMD simulation of central nuclear 

collisions, ^,-(g+u) lo<s'<i/ m~{s+v) l^- - ^■^^■ 

E. The V' to J/tl; ratio 



Figures |T0| and |TT] show the ^jj' /J/ip ratio in pA and AB collisions as a function of 
AB. UrQMD calculations in Scenarios II and III are compared with the experimental data 
The observed weak A dependence of the ratio in pA reactions is consistent with the 



hypothesis of the formation of a 'preresonance' state, e.g. |jTI|]. If the charmonium states 



are initially produced in the \ccg) configuration and color neutralization takes sufficiently 
long, the produced J/ip^s and ip'^^ are equally suppressed by nuclear absorption. This is 
reflected in Fig. |TD|, showing the calculations within Scenario II. In pA collisions, comover 
absorption is only a minor contribution to the total charmonium absorption. Consequently, 
the universal nuclear absorption cross section of the color-octet states results in a constant 
ifj' / J /ip ratio. 

In Scenario III, the absorption cross sections for the individual charmonium states evolve 
linearly with time. Thus, the ip' , Xcio and Xcii states have absorption cross sections similar 
to that of the J/il) as long as r* < tJ . Only for r* > tJ the cross sections differ strongly, 
because then a{J/ipN) has reached its asymptotic value while the cross sections of the larger 
resonances still increase. The differences in the calculated J/ip and ip' survival probabilities 
for pA reactions, as shown in Fig. |TT], are therefore mainly due to nuclear dissociation 
processes in the later stage of the reaction and interactions with secondaries. Even in pA 
reactions comover interactions are not negligible in Scenario III because leading mesons are 
taken into account. 



Figures |T0| and |r^ also show the calculated ip' /J/ijj ratio in S+U and Pb+Pb reactions. 
In Scenario II, the ijj' /J/ip ratios in S+U and Pb+Pb reactions very roughly agree with the 
experimental values. In Scenario III, dissociation by mesons leads to an even stronger de- 
crease of the %l)' / J/tl) ratio in S+U and Pb+Pb collisions than seen in the data 0. Although 
shortly after charmonium production the charmonium- comover dissociation cross sections 
are similar for all states, the kinematic threshold leads to different dissociation probabilities. 

The ijj' / J /ip ratios in nucleus- nucleus collisions seem to favor Scenario II as compared 
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to Scenario III, although neither describes the data satisfactorily. However, the parameter 
space of the models has certainly not yet been fully explored. Moreover, for a complete 
understanding of the ratio, quantum interference effects |]5^ as well as refeeding processes 



such as 'nJ/il) -^ ip'iT [^|58| must be considered. As can be seen in Fig. [TT|, the ip' survival 
probability in S+U collisions in Scenario III would need to be increased by more than 50% 
in order to agree with the experimental observations. To test the effect of refeeding, we 
have implemented the process ttJ/ijj — > i/j'tt in UrQMD with an energy independent cross 
section above the kinematic threshold, ii^th = ^^tp' + "^tt- The refeeding cross section is 
assumed to have a time dependence like that of the nuclear and comover absorption cross 
sections. In this case, a.„ji^^^i.„ = 7 mb would increase the ip'/J/ip ratio by 50%. This 
cross section seems unrealistically large because this exclusive channel a is larger than the 
inclusive cr(7r J/ip) = 2.4 mb we have used elsewhere even though the mass threshold is larger 
than 71 J/ip -^ DD. In order to further explore the effectiveness of possible ip' refeeding 
processes, we have implemented the channel MJ/ip — > ip'M, where M stands for any meson 
species. Then the refeeding cross section must be ctmj/ip^iP'm = 3 mb in order to increase 
the V' yield by 50%. 

Finally, note that a change of the ip' dissociation cross section alone evidently has a 
direct effect on the ip'/J/ip ratio while the J/ip survival probability itself is hardly affected. 
Therefore, in view of the theoretical uncertainties concerning charmonium dissociation cross 
sections and formation times it seems premature to draw definite conclusions from the model 
comparison presented here. 



F. Pb+Pb at 160 GeV vs. 80 GeV 

Since lowering the SPS beam energy has been discussed, it is worthwhile to check if such 
a change could distinguish between an evolving cc dissociation cross sections and a constant 
cross section. With a constant cross section, nuclear absorption should be insensitive to 
the 7 factor, i.e. the Lorentz contraction, of the nuclei because the reduced path length is 
canceled by the enhanced density. There is no cancelation, however, in the case of evolving 
cross sections. In the limit of very high beam energies the nuclei become transparent for 
charmonium states despite the large nuclear density, since the absorption cross section is es- 
sentially zero immediately after production. Comover absorption, on the other hand, should 
increase with energy in both scenarios because the number of secondaries rises. However, the 
average comover velocity should be also enhanced which increases their average formation 
times in the center of mass frame counteracting the effect of increasing density. Thus for 
very high energies, disregarding leading hadrons, absorption by comovers must eventually 
become negligible due to the large formation times. In any case, the comover density and its 
time evolution varies with energy and can lead to different characteristic behaviours within 
the two scenarios. 

We have compared the calculations with the standard choice of parameters for the calcu- 
lation with constant cross sections. Scenario II, and the evolving cross sections. Scenario HI. 
The resulting survival probabilities for minimum bias Pb-|-Pb collisions at E'lab = 160 GeV 
and 80 GeV are presented in Table |. In fact, the calculated energy dependence of the nu- 
clear absorption within the two scenarios behaves as anticipated. Dissociation by nucleons 
stays constant in Scenario II while it decreases with increasing beam energy in Scenario HI. 
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The numerical results show that the comover absorption gets slightly stronger with increas- 
ing energy both for constant and evolving dissociation cross sections. The combined effects 
of absorption by baryons and mesons lead to a decrease of the J /ip survival probability for 
80 GeV as compared to 160 GeV in the case of evolving cross sections and an increase for 
constant cross sections. However, the relative change of the survival probability is only on 
the order of a few percent in both cases. Thus the experimental distinction between the two 
scenarios with the suggested change of beam energy is unlikely. 

G. Rapidity dependence of J /ip suppression 

Another means by which the charmonium production mechanism might be determined is 
the study of the survival probability as a function of rapidity. We investigate pA interactions 
because absorption effects are limited to only one nucleus and comover suppression is signif- 
icantly reduced in pA relative to AB collisions. A constant absorption cross section should 
be independent of rapidity while an evolving cross section should produce more absorption 
close to the target since the asymptotic state is more likely to interact inside the target. 



Figure |T^ shows the rapidity dependence of the J/ip survival probability in p+U reactions 
at -Eiab = 200 GeV. The predictions of Scenario II and Scenario III are compared. The 
charmonium absorption is dominated by nuclear absorption in this system although comover 
interactions are not completely negligible, in particular in Scenario III, see Figs. and |^. 
The scenario of a constant color-octet dissociation cross section leads to a constant rapidity 
dependence of the J/tp survival probability because the effective path length of charmonium 
'preresonance' states does not depend on their longitudinal velocity as long as yx{cc) > 1/target- 
The assumption of evolving color-singlet states leads to an increasing survival probability 
with increasing rapidity. At forward rapidities, the larger 7 factor of the charmonium states 
leads to an increased formation time in the NN center of mass frame. The evolution of the 
dissociation cross sections is therefore delayed. The charmonium states suffer effectively less 
absorption during the passage through the target nucleus. The reverse effect is limited by the 
asymptotic values of the absorption cross sections. Therefore, the J/ip survival probability 
saturates at negative rapidities. 

The experimental examination of the rapidity dependence in pA reactions is a very 
promising test of the production mechanism because of the very distinct behaviour of the 
J/ip survival probability. It would be particularly interesting to carry out this experiment 
at about the same energy as the nuclear collisions studied by NA38 and NA50. Detailed 
analyses of existing experimental data on the xp dependence of J/tp production in pA 
reactions [P^,P^ should also be used to test the vahdity of different production mechanisms 



as done in Ref. |5^. There, it is also shown that the shape of the xp dependence changes with 
energy. It is conceivable that a realistic scenario of charmonium production and absorption 
involves a superposition of color-singlet and color-octet contributions ^^, however, the 
relative importance may well be energy dependent. 
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H. A closer look into interactions with secondaries 



From the results presented in Sec. |111 C| we conclude that the observed J /ip production 



cross sections in pA and AB collisions can be explained in purely hadronic scenarios provided 
that secondary particles can effectively dissociate charmonium. In Scenarios II and III we 
have employed universal and energy independent^ comover absorption cross sections for the 



individual cc meson states, see Fig. |1^. Given the X{cc)N cross sections ||23[, the size 
of the comover cross sections — with a reduction of F = 1/2 and 2/3, respectively — 
are presumably of the right order, if similar collision energies are considered. However, 



calculations within the framework of short-distance QCD [Q show that J/ip dissociation 
by hard gluons is strongly damped close to the kinematic threshold. Such a calculation 
is shown in Fig. O for the process J/^tt. This would significantly increase the survival 



probability. However, the comover absorption would not be completely suppressed because 
of the feeding from the less bound Xc and ip' states. For example, the survival probability due 
to interactions with comovers in minimium bias Pb(160 GeV)+Pb collisions in Scenario II 
increases from Sm = 0.61 to 5*^ = 0.69 if the J/ip itself cannot be dissociated by mesons. 
As a result, the total survival probability in this system would increase from Stot = 0.31 to 
^tot = 0.35. 

Calculations of J/ip it and J/ip p dissociation cross sections in the framework of a meson 
exchange model also show a suppression at low energies even for the exothermic J/ip p -^ DD 
process |^^. The cross section for the dissociation by pions according to these calculations 



is shown in Fig. |T3|. Note that the inclusion of a form factor would further decrease the cross 
sections [p4| . 

From phase space arguments alone, i.e. assuming an energy independent matrix element, 
one would infer an enhancement for low relative velocities in the incoming channel if the 
dissociation process is exothermic. Calculations of cross sections involving heavier mesons 
in the incoming channels within meson exchange models would render valuable information 
in order to estimate comover absorption probabilities. Also, a deeper understanding of the 
reactions involving the Xc and ip' states is needed. 

In a chemically equilibrated environment at temperatures less than 300 MeV, the con- 
sideration of interactions with vr's and p's might be sufficient to quantitatively estimate the 
total comover absorption. In any kind of thermal hadronic scenario, higher resonances are 
negligible due to the low density, pu ~ exp(— m/^/T). However, Table |I| shows the rela- 
tive importance of the mesons which contribute most to J /ip dissociation according to the 
nonequilibrium UrQMD calculation. Scatterings of the J/^/^'s with vr's and p's are indeed 
the dominant dissociation processes. However, together they are responsible for only about 
40% of the total comover absorption. The composition is similar in Scenarios II and HI. 
Twenty channels are left out of Table |I|, each of them contributing ?« 1% or less. However, 
they account for about 15% of the total absorption. 

The string fragmentation scheme employed in UrQMD [Q determines the relative 



abundancies of hadron states in the nonequilibrium simulation. The choice of parame- 
ters is far from unique since the model can be compared to data only for experimen- 



^ Above the respective kinematic dissociation threshold, i.e. (Tdiss(\/s) = S{y/s — £'th)o"o 
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tally identifiable hadrons. In UrQMD, all states of eight different meson multiplets 
(J^*" = 0""*", 0+"*", 1 , l+~, 1++, 2++ plus two excited 1 multiplets) may be populated 
in a string fragmentation process. The probability to form a meson from one of these mul- 
tiplets is chosen to be proportional to the spin degeneracy and inversely proportional to the 
average mass of the multiplet |^ . 



In Fig. [Tj, the average J /ip n collision energies as functions of time are shown for three 
different scenarios: constant cross sections including interactions of leading mesons (c/. Ref. 
||60|), constant cross sections excluding leading mesons (Scenario II), and evolving cross 
sections sections including interactions of leading mesons (Scenario III). We studied central 
Pb(160 GeV)+Pb collisions. All model simulations yield high average collision energies of 
(E) ^ 4.8 GeV at t ^ 1 fm/c. As stated in Ref. [^, these collision energies are not 



compatible with an equilibrated hadron gas scenario since they correspond to a meson gas 
with a temperature oiT ^ 750 MeV. The collision energies and the scattering rates at later 
times, t > 2 fm/c, are both consistent with temperatures of T = 140 ± 20 MeV [^], thus 
allowing for the application of a thermal reaction kinetics model. 

An important difference between the three scenarios presented in Fig. |T^ is shown in 
Fig. ^ which determines the impact of the highly nonthermal secondaries at early times on 
the total comover absorption. The scattering rate as a function of time, normalized to one 
J / ip-Taeson collision, is depicted. It can be seen that at early times t < 2 fm/c, dissocia- 
tion by comovers is almost exclusively due to interactions of leading mesons. Moreover, if 
leading meson interactions are taken into account, they clearly dominate the total comover 
absorption. This is true for the scenario with constant as well as with evolving cross sec- 
tions. Note that the dissociation cross section of a leading meson is only 1/2 of that of 



a fully formed meson. As can be seen from Figs. |T^ and |T^, the applicability of thermal 
momentum distributions for the comover-charmonium interactions in heavy ion collisions 
depends heavily on the underlying model assumptions. A thermal comover model seems 
to be a reasonable approximation in Scenario II while it cannot be justified in Scenario III 
where leading mesons are included. 

Figure |16| shows the time integrated collision spectrum of J/ip vr interactions in central 
Pb(160 GeV)+Pb reactions according to the UrQMD simulations. The high energy tails of 
the spectra calculated within the different model scenarios reflect the interactions at early 
times. They correspond to temperatures of T ^ 300 MeV. However, it is also clear from 
Fig. 0, that for all presented scenarios, more than 75% of J/ip vr interactions occur at 
collision energies below 4.25 GeV. Thus, the threshold behaviour of all involved dissociation 
channels will finally decide, whether comover absorption is effective enough to explain the 
suppression of charmonium yields in nuclear reactions. 

IV. CONCLUSIONS 

We have studied the production and absorption of charmonium states at SPS ener- 
gies using a microscopic hadronic model, UrQMD. The model simultaneously reproduces 
the final rapidity distributions of produced hadrons and net-baryons in pp and AB colli- 
sions. The simulations also agree with the measured total transverse energy spectrum for 
S(200 GeV)+Au (NA35) and Pb(160 GeV)-fPb (NA49) as well as with the Et dependence 
of the Drell-Yan muon pairs in Pb(160 GeV)+Pb (NA50). At the same time the calcu- 
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lations are in significant disagreement with the S(200 GeV)+U Et dependence of dimuon 
production (NA38). 

Our study shows that nuclear absorption alone cannot yield prefect agreement with all 
measured J/ip cross sections in pA and AB collisions at SPS energies. This is in accord 
with the findings of previous studies with semianalytical models. However, we find that a 
global fit to the AB systematics with a nuclear absorption cross section of 9 mb is not yet 
completely ruled out by the data taking into account its systematic and statistical errors 
and its in part still prelimary status. 

The UrQMD model renders better agreement with the J/ip data if comover dissociation 
is included and the charmonium is produced in color-singlet states which evolve in time. In 
this case, all parameters were fixed by other data before comparison to the NA50 data. The 
Et dependence of the J/'j/'/Drell-Yan ratio is also reproduced in this model. The ip'/J/ijj 
ratio agrees with the pA data, while the ip' suppression in S+U and Pb+Pb is overestimated. 
Assuming produced color-octet states with a time independent, universal cross section and 
comover absorption of the fully formed cc states, yields equally good agreement with the 
J/ip data. The overall trend of the ip' /J/ip systematic for different projectile and target 
combinations is fairly well reproduced, although significant quantitative deviations from the 
data are found in nuclear collisions. 

The rapidity dependence of the J/ip survival probability in pA reactions is a promis- 
ing observable by which the contributions of color-singlet and color-octet states to nuclear 
absorption processes can be disentangled. 

The effectiveness of charmonium dissociation by interactions with secondaries is appar- 
ently crucial for the understanding of the experimental data in hadronic scenarios. The 
microscopic simulation shows that the charmonium-comover interactions are not close to 
an equilibrium situation. Heavy meson resonances contribute to a large fraction of all 
charmonium-comover interactions. The calculation of energy dependent cross sections for a 
large number of possible hadronic dissociation channels is a prerequisite for a quantitative 
understanding of comover absorption. 
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TABLES 

TABLE L Comparison of J /ip survival probabilities for minimum bias Pb+Pb collisions at 
different energies for constant and evolving dissociation cross sections, i.e. Scenario II and Sce- 
nario III, respectively. 



-E'lab 


all 


all 

'JJVI 


all 

'^tot 


^B 


qlll 


qlll 

'-'tot 


160 GeV 

80GeV 


0.51 
0.51 


0.61 
0.64 


0.31 
0.33 


0.57 
0.50 


0.49 
0.51 


0.33 
0.32 



TABLE II. Contributions of different meson species to the total J/tp comover absorption in 
central Pb+Pb collisions at 160 GeV according to the UrQMD simulation with Scenario II and III. 
The percentage of comover absorption due to the most dominant channels are shown. 



Species 


Scenario II 


Scenario III 


vr 


19.7 ±0.3 


18.6 ±0.3 


P 


20.8 ± 0.4 


18.6 ±0.3 


K 


11.6 ±0.3 


11.1 ±0.3 


K*{892) 


8.9 ±0.3 


8.1 ±0.2 


7] 


5.4 ±0.2 


6.4 ±0.2 


LO 


8.4 ±0.3 


6.3 ±0.2 


02(1320) 


3.1 ±0.1 


4.4 ±0.2 


ai(1260) 


2.6 ±0.1 


4.2 ±0.2 


6i(1235) 


3.6 ±0.2 


3.8 ±0.2 


ao(980) 


3.0 ±0.1 


2.9 ±0.1 


Sum: 


87.1 ±0.8 


84.4 ± 0.7 
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FIG. 1. Differential cross section of the transverse energy produced in Pb(160 GeV)+Pb and 
S(200 GeV)+Au collisions within the pseudorapidity range 2.1 < rj < 3.4. The UrQMD results are 



compared with data from |42|. 
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FIG. 2. a) Number of Drell-Yan muon pairs in Pb+Pb as a function of the produced neutral 
transverse energy within 1.1 < rj < 2.3 with arbitrary normalization. The UrQMD result is shown 
with data from NA50 E^ where the S-p-axis of the data is rescaled by 0.8. The modification 
is motivated by the recently published Et — Ezdc contour plot Q compared to the previously 
published analysis [^Q. b) Number of Drell-Yan muon pairs in S+U as a function of the produced 
neutral transverse energy within 1.7 < rj < 4.1 with arbitrary normalization. The UrQMD result 



and data from NA38 [44| are shown. 
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a) 



b) 
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FIG. 3. Time evolution of a central Pb+Pb collision at -Eiab = 160 GeV according to the 
UrQMD simulation with time independent absorption cross sections. The interaction points of 
charmonium states in the z-t-plane are shown, a) The cc production points, b) The cc dissociation 
points due to interactions with nucleons, if the momentum distribution of the charmonium states 
and the nuclear stopping is neglected, c) The momentum distribution of the charmonium states is 
now included while the nuclear stopping is still neglected, d) The momentum distribution of the 
charmonium states and nuclear stopping is taken into account. 
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FIG. 4. Time evolution of a central Pb+Pb collision at £^iab = 160 GeV according to the 
UrQMD simulation with time independent absorption cross sections. The interaction points of 
charmonium states in the z-t-plane are shown, a) The cc dissociation points due to interactions 
with mesons including leading hadrons. b) The cc dissociation points due to interactions with 
mesons without leading hadrons. 
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FIG. 5. J/V'-production cross sections times dimuon branching ratio in the kinematical domain 
< Vcm < 1 and | coshes I < 0.5, and rescaled, if necessary, to piab = 200 GeV as a function of 



AB. Open circles denote the semi-analytical Glauber calculations of |11| while the results of the 
microscopic UrQMD simulation without nuclear stopping and momentum distribution of the J/ip, 
Scenario I, are shown as full circles. A universal and constant absorption cross section of a = 7.3 mb 
is used for both. The data (open triangles) are from O. 
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FIG. 6. J/V'-production cross sections times dimuon branching ratio in the kinematical domain 
< Vcm < 1 and | cos^csl < 0.5, and rescaled, if necessary, to piab = 200 GeV as a function of 
AB. The results of the microscopic UrQMD simulation without nuclear stopping and momentum 
distribution of the J/il^, Scenario I, are shown as full circles. A universal and constant absorption 
cross section of o" = 9 mb is used. The J/tp production cross section in p+p reactions (the absolute 
normalization) is arbitrarily chosen. The data (open triangles) are from H. 
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FIG. 7. J/V'-production cross sections times diinuon branching ratio in the kinematical domain 
< Vcm < 1 and | cos 6cs\ < 0.5, and rescaled, if necessary, to piab = 200 GeV as a function of AB. 
Full circles denote the results according to Scenario II. Nuclear absorption is calculated without 
stopping using a universal dissociation cross section of cr = 4.8 mb, cf. Ref [^. Comover absorption 
is calculated in a full UrQMD simulation without leading hadrons. The charmonium-meson cross 
sections {X{cc) + vr, X{cc) + p, etc.) are reduced by a factor of 1/2 from the cross sections 
for X(cc) + N from Ref. pS]. Open circles denote the production cross sections if only nuclear 
absorption is considered. The data (open triangles) are from ||3|. 
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FIG. 8. J/^-production cross sections times dimuon branching ratio in the kinematical domain 
< Vcm < 1 and | cos 6cs\ < 0.5, and rescaled, if necessary, to pi^h = 200 GeV as a function of AB. 
The UrQMD calculation with evolving dissociation cross sections, Scenario III, is shown. Open 
circles denote the production cross sections if only nuclear absorption is considered. The data 
(open triangles) are from Q . 
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FIG. 9. The ratio of J /ijj to Drell-Yan production as a function of Et for Pb+Pb collisions 
at 158 GeV (a) and S+U at 200 GeV (b) according to Scenario III. The experimental data are 
from Ref. [^ and [^, respectively. The normalization factor Ba^/a^^ = 46 in pp interactions at 
200 GeV is taken from Ref. |11|. An additional factor of 1.25 |]l5[ has been applied to the Pb+Pb 
calculation in order to account for the lower energy (the J /ij: and the Drell-Yan cross section have 
a different energy dependence). The Et values of the experimental S+U data from [^ are scaled 
by a factor 1.25 assuming that the UrQMD calculation reflects the correct Et distribution, see 
Fig. I 
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AB 



FIG. 10. ip' / J /ip ratio in pA, S+U and Pb+Pb interactions as function of AB for universal 
nuclear dissociation cross sections and distinct charmonium states interacting with comovers after 
their formation time, Scenario II. The data (open triangles) are from |^,5(:]. 
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FIG. 11. 'ip' / J /tp ratio in pA, S+U and Pb+Pb interactions as function of AB for evolving 



dissociation cross sections, Scenario III. The data (open triangles) are from |^,5(:]. 
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JZ-^ suppression in p(200 GeV)+U 
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FIG. 12. J/f/'-survival probability in p(200 GeV)+U reactions as a function of rapidity. Calcu- 
lations with Scenario II (open circles) and Scenario III (full circles) are shown. 
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short-distance QCD 
meson exchange 




10 



FIG. 13. J I'^'K dissociation cross section as calculated in short-distance QCD according to |34] 
in comparison with a calculation with a meson exchange model [^^. Also shown are the values 
of the J/'0-meson absorption cross section chosen for the model calculations with evolving cross 
sections, a = 2.4 mb (Scenario III), and for the calculations with constant cross sections, color-octet 
nuclear absorption and no leading mesons, a = 1.8 mb (Scenario II). 
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FIG. 14. Average J/tpir collision energies in central Pb(160 GeV)+Pb reactions as a function 
of time for different comover scenarios. We show the distributions for constant dissociation cross 
sections including leading meson interactions (full line), excluding leading mesons, Scenario II 
(dashed line) and for evolving cross sections with leading mesons, Scenario III (dotted line). 
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FIG. 15. J/^-meson collision rates, normalized to one dissociation, in central Pb(160 GeV)+Pb 
reactions as a function of time. We show the results for constant dissociation cross sections including 
leading meson interactions (full line), excluding leading mesons. Scenario II (dashed line) and for 
evolving cross sections with leading mesons, Scenario III (dotted line). 
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FIG. 16. Time integrated J/tpir collision spectrum in central Pb(160 GeV)+Pb reactions ac- 
cording to the UrQMD calculation (normalized to one collision) in comparison to a thermal hadron 
gas. Shown are the energy spectra resulting from constant dissociation cross sections including 
leading meson interactions (full line), excluding leading mesons, Scenario II (dashed line) and 
for evolving cross sections with leading mesons. Scenario III (dotted line). The J/iJ)tt collision 
spectrum in an ideal gas of tt's at T = 300 MeV is indicated by the dash-dotted curve. 
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